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The ground-state properties of the t — J model on a d-dimensional hypercubic lattice are examined 
in the limit of large d. It is found that the undoped system is an ordered antiferromagnet, and 
that the doped system phase separates into a hole-free antiferromagnetic phase and a hole-rich 
phase. The latter is electron free if J > it and is weakly metallic (and typically superconducting) if 
J < At. The resulting phase diagram is qualitatively similar to the one previously derived for d = 2 
by a combination of analytic and numerical methods. Domain wall (or stripe) phases form in the 
presence of weak Coulomb interactions, with periodicity determined by the hole concentration and 
the relative strength of the exchange and Coulomb interactions. These phases reflect the properties 
of the hole-rich phase in the absence of Coulomb interactions, and, depending on the value of J/t, 
may be either insulating or metallic (i.e. an "electron smectic"). 



In this paper, the zero-temperature properties of 
the t — J model of a doped antiferromagnet on a d- 
dimensional hypercubic lattice are evaluated using a sys- 
tematic expansion in powers of 1/d. For each property of 
interest the leading behavior in the large d limit is com- 
puted, and in some cases, just to prove how tough we are, 
corrections up to order 1/d 5 are obtained. These results 
are obtained by breaking the full Hamiltonian into an un- 
perturbed piece, Hq, and a perturbation, Hi, and then 
reorganizing conventional perturbation theory in powers 
of Hi into a 1/d expansion. Of course the partition of the 
Hamiltonian may be chosen for calculational convenience, 
since it does not affect the results. The convergence of 
this expansion will not be addressed, although we believe 
it to be only asymptotic. 

Our procedure differs from the extensive recent work 
on the related problem of the Hubbard and Falicov- 
Kimball models in large dimensionEI in the way the large 
dimension limit is taken. First of all, we do not assume 
that the ratio J/t of the exchange integral J and the 
hopping amplitude t is parametrically small as d — > oo. 
The previous studies assumed that t is proportional to 
1/Vd so that, when J is expressed in terms of the onsite 
interaction U, it follows that J/t = At/U ~ 1/Vd. (The 
phase diagram will be studied for parametrically small 
values of J/t in Sec. VIII, but our results are less com- 
plete in this case, because of the difficulty of controlling 
perturbation theory in this limit.) Secondly, the hyper- 
cubic lattice is bipartite, i.e. it can be broken into two 
sublattices, which we label "black" and "red" , such that 
the Hamiltonian has interactions only between sites on 
different sublattices. This favors the classical Neel state, 
which has a uniaxial magnetization with opposite sign on 
the two sublattices. By contrast, earlier studies, which 
were primarily concerned with the Mott transition and 
possible non-Fermi liquid states of the Hubbard model, 
assumed a non-bipartite lattice which frustrates the Neel 



state. For both reasons, this previous work does not shed 
much light on the behavior of doped antifcrepmagnets. A 
notable exception is the work of van Dongeno on the small 
U limit of the Hubbard model on a hypercubic lattice 
which found, as we do, that the weakly-doped antiferro- 
magnetic phase is unstable to phase separation, even if 
the parameters are scaled as J eff ft = At/U ~ 1/Vd. 

Throughout this paper, units are chosen such that the 
lattice constant, h, and Boltzmann's constant are all 
equal to one. 



I. SUMMARY OF RESULTS 
A. Results in Large Dimension 

Our principal result is the global zero-temperature 
phase diagram as a function of J/t and hole concentra- 
tion x, in the limit of large d, as shown in Fig. 1. It 
is immediately clear that in most of the phase diagram, 
the undoped (ordered) antiferromagnetic phase coexists 
with a hole-rich phase. For J/t > 4, the hole-rich phase 
is electron free; otherwise it contains an exponentially 
small but non-vanishing concentration of electrons. In 
the intermediate coupling regime, 2 < J/t < 4, the resid- 
ual attraction ( J ) between electrons is great enough to 
overcome the hard-core repulsion, and leads to a BCS 
instability of the dilute metal, producing an s-wave su- 
perconducting state at exponentially low energy scales. 
At smaller values of J/t, the net interaction between 
electrons is repulsive. This implies that the system ei- 
ther remains metallic down to zero temperature or ex- 
hibits higher-angular-momentum pairingtj via the Kohn- 
Luttinger mechanism. □ 

A peculiarity of the phase diagram in Fig. 1 is that the 
boundary of the two-phase region intersects the J/t = 
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FIG. 1. Phase diagram of the t — J model in the limit 
d — > oo: Here x is the hole concentration (1 — x is the electron 
concentration). The phase boundary is given by Eq. ([53|), 
artificially setting d = 2. "Two-phase" labels the two-phase 
region, where a uniform density phase is thermodynamically 
unstable, "SC" labels a region of s— wave superconductivity, 
and "M" labels a region of metallic behavior with repulsive 
interactions, which presumably has an ultra-low temperature, 
superconducting instability due to the Kohn-Luttinger effect.u 

axis at a non-zero value of x. This is not likely to be 
correct in any finite dimension. For small x and large but 
finite dimension, we expect that in the limit J/t — > 0, the 
ground state is a ferromagnetic Fermi liquid, and hence 
the model does not phase separate. In Section VIII, wc 
discuss the behavior of the model for J/t parametrically 
small, J/t ~ 1/y/d. Here the l/d expansion is slightly 
more difficult to control, so our results, summarized in 
Fig. 2, are incomplete. The resulting conjectural phase 
diagram for large but finite d embodies all the insights 
gained from studying the d — > oo limit, but corrects the 
unphysical features of the phase diagram in Fig. 1. 

We have also studied the behavior of one or two doped 
holes and the character of charged domain walls in the 
antiferromagnet. It will be seen that the latter are sta- 
bilized by a long-range Coulomb interaction. These 
studies bring out an important characteristic of our large- 
d expansion. Whenever a hole moves in the antiferromag- 
netic background, it may break a number of bonds of or- 
der d at each hop. Consequently, for such processes, the 
physics is exchange dominated for large-<i and it amounts 
to an expansion in powers of t/ J. This is true of the mo- 
tion of one or two holes and of domain wall fluctuations in 
which holes hop into the environment. However the ques- 
tions of phase separation, domain wall phase equilibrium, 
and superconductivity at low electron concentration are 
not subject to this limitation. 

The ensuing discussion will be organized by order of 
increasing hole concentration. 

To leading order in 1 / d the states of minimum energy 
of a single hole lie precisely on the magnetic Brillouin 
zone which also is the Fermi surface of the noninteract- 
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FIG. 2. Conjectured Phase diagram of the t — J model for 
large but finite d: This figure should be viewed as a blowup of 
the small J/t portion of Figure 1. The horizontal line repre- 
sents the small J/t extension of the phase boundary in Fig. 1; 
in fact, in large d, this line would be exponentially close to the 
top of the figure, but we have drawn it, as in Fig. 1, at a posi- 
tion obtained by setting d — 2 in the large d expression. The 
boundary of the fully polarized ferromagnetic metallic phase 
(labelled "F") is drawn in accord with the large d expression 
in Eq. (frH| There might be other, lower energy phases, (e.g. 
high-density stripe phases) that could occur below these two 
phase boundaries, in the region marked "two-phase", espe- 
cially close to the point of intersection. 



ing system with a half-filled band. The massive degen- 
eracy of these low energy hole states is lifted by terms 
of 0(l/d 4 ), and it is found that the absolute minimum 
occurs at k = (tt/2) < 1, 1, 1, ... > together with points 
related by the point-group symmetry. Moreover, as de- 
duced previously by TrugmanQ in studies of two holes in 
a two-dimensional antiferromagnet, we find that propa- 
gation of pairs of holes is no less frustrated than is the 
propagation of a single hole, because of a subtle effect 
of Fermi statistics. There is, however, an effective at- 
traction ~ l/d between two holes due to the fact that 
two nearest-neighbor holes break one less antiferromag- 
netic bond than two far-separated holes; this attraction 
always leads to a two-hole bound state. 

An interesting metastable state is a charged magnetic 
domain wall (i.e. a d — 1 dimensional hypersurface with 
finite hole concentration and suppressed magnetic order) . 
We have found that the most stable domain wall has an 
electron density which is, to leading order in l/d, equal 
to that of the hole-rich phase which can exist in equilib- 
rium with the antiferromagnet. Thus domain walls can 
be viewed as a form of local phase separation. Also the 
domain wall configuration with the lowest surface ten- 
sion (i.e. energy per unit hyperarea of wall) is the "ver- 
tical" site-centered 7r (antiphase) discommensuration in 
the antiferromagnetic order; i.e. it is parallel to a single 
nearest-neighbor vector and odd under reflection through 
a site-centered vertical hyperplane. 
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We have considered the effect of weak, long-range 
Coulomb interactions as a perturbation. While this 
study is not exhaustive, we conclude that, for a sub- 
stantial range of parameters, the ground state consists 
of a periodically-ordered array of optimal domain walls 
of the sort described above, especially when x is small 
but not too small. In this range of x, the ground state is 
insulating for J/t > 4, and metallic, for J/t < 4. The lat- 
ter phase is an "electron-smectic"Ll which exhibits crys- 
talline order in one direction and liquid-like behavior in 
the transverse (d— 1) directions. The liquid features are 
associated primarily with the motion of electrons along 
the domain wall, and they may be metallic or condensed 
into a superconducting state. 

We have argued previously that the competition be- 
tween a local tendency to phase separation in a doped 
antiferromagnet and the long-range Coulomb repulsion 
between holes produces a large variety of intermediate 
scale structures, including arrays of domain walls, which 
are significant features of doped antiferromasaiets that we 
have called "frustrated phase separation" Wu. However, 
these phenomena have not previously been derived from 
a microscopic magnetic model.El It is particularly striking 
that, in the appropriate range of parameters, charge and 
spin density wave order coexist with metallic, and even 
superconducting behavior. 

B. How Large Are d = 2 and d — 3? 

Large d is, of course, only of academic interest; we 
are interested in the physical dimensions, d = 1, 2, and 
3. The properties of the oiteydimensional electron gas 
(1DEG) are well understoocfcJ by now, and exhibit be- 
havior that is quite dimension specific. Moreover, for 
most of the conceivable ordered states, the lower critical 
dimension for long range order at zero temperature is one, 
so the 1DEG is not likely to be well understood in terms 
of adiabatic continuity from large dimension. However, 
long range order at zero temperature is quite robust in 
both two and three dimensions, so there is every reason 
to expect that a l/d expansion will capture the essential 
physics of many of the zero-temperature thermodynamic 
states. 

To test this conjecture, we would like to make both 
qualitative and quantitative comparisons between the re- 
sults of the large d theory and any available exact, or 
well-controlled numerical or analytic results in two and 
three dimensions. Table 1 gives a quantitative compari- 
son between the 1/d expansion and well-established nu- 
merical results for the undoped system, i.e. for the spin- 
1/2 Hcisenberg antiferromagnet. It can be seen that the 
ground-state energy can be obtained from the low-order 
expansion in powers of 1/d to 0.6% accuracy or better. 
By carrying the series to higher order, and possibly do- 
ing a Pade analysis of the series much improved accuracy 
for all physical quantities could be expected. In Sec. XI 



comparison will be made between numerical results and 
the results of perturbation theory about the Ising limit 
(Table 2). 
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Table 1: .Comparison of the results of exact numer- 
ical studiesEiJ (the row labeled "exact") on the two- 
dimensional spin-1/2 Heisenberg antiferromagnet with 
the perturbative results in powers of 1/d derived in the 
present paper. (We have been unable to find correspond- 
ing "exact" three dimensional results.) The dimension is 
indicated by the arguments of the computed quantities. 
The rows labeled "Upper" and "Lower" give the rigorous 
upper and lower bounds on the energies obtained in the 
text. The approximate results are obtained by setting 
y = J±/J z = 1, V = 0, and d = 2 or 3 in the series ex- 
pansion, evaluated to the stated order. All energies are 
measured in units of J/d, and the magnetization m is 
quoted in units in which gfig = 1, where hb is the Bohr 
magneton. 
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FIG. 3. Zero Temperature Phase Diagram of the two di- 
mensional t — J model, deduced from numerical studies of 
finite size systems with up to 60 electrons, as well as from 
various analytic result** This figure is abstracted from Hell- 
berg and ManousakisJij 

Qualitative comparisons can be made with the phase 
diagram of the two dimensional t — J model which has, 
been deduced from combined analytic and numericallijEll 
studies. Figure 3, abstracted from the work of Hell- 
berg and ManousakisE-2 shows the phase diagram deduced 
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from numerical studies of systems with up to 60 electrons. 
As in large d, there is no thermodynamically stable zero- 
temperature phase with dilute holes for any J / t. Indeed, 
aside from the behavior of the boundary of the two-phase 
region at very small J/t, the phase diagrams in Figs. 1 
and 3 are similar. As suggested above, when the patholo- 
gies of the formal d — > oo limit are removed by taking into 
account the new processes that become important at pa- 
rametetrically small values of J/t ~ d^ 1 / 2 , one obtains 
for large but finite d the phase diagram shown in Fig- 
ure 2, which is topologically equivalent to Figure 3. (Of 
course, in d = 2, parametrically small values of J/t are 
not all that small, so there is no reason to expect quan- 
titative agreement with the large d results. The critical 
value of J/t = Y c at which the phase-coexistence line 
deviates from x = 1 ist-3 Y c — 3.4367 in d = 2, and is 
rather well approximated by the value Y c — > 4 as d — > oo. 
However, the slope of the phase coexistence line in d = 2 
is much steeper than would be deduced from the large d 
theory.) Similar detailed information on the three dimen- 
sional t — J model is not available at this time, although 
arguments presented hithertoOjij suggest that the phase 
diagram is qualitatively similar to that in d = 2, consis- 
tent with the expectations from the 1/d expansion. 

Our calculation of the spectrum of one hole in an an- 
tiferromagnet may be compared to the numerical calcu- 
lations of Dagotto et alx3 on d = 2 systems with 16x16 
sites and J/t = 0.4. They found that the one-hole spec- 
trum is well represented by the two-dimensional version 
of the expression in Eq. ( |3q ) , confirming the qualitative 
accuracy of the large d expression. However the values 
of the parameters obtained to leading order in 1/d are 
quantitatively quite far from the exact results, and this 
discrepancy is made worse by the inclusion of higher or- 
der terms. This is not unexpected in view of the fact 
that large d drives the motion of a single hole into the 
exchange-dominated limit. In particular, it is clear from 
Eq. (B7|) that the large-<i expansion gives a negative 
value for the bandwidth W in d = 2, unless J > 0.93t. 
Thus it is essential to compare the large-d expansion to 
numerical results at large J/t. Specifically, from Eq. 
( |36| ) with y = 1, the bandwidth for d = 2 is given by 
dW/t = 2.125t/J-l.83(t/J) 3 . It would be interesting to 
compare this result with numerical calculations for large 
J/t extrapolated to the thermodynamic limit. Martinez 
and HorschEJj have found that an approximate treatment 
of the motion of a single hole gives dW/t = 2t/J for large 
J/t, which agrees very well with our large-<i result!— Via a 
variational calculation, Bonisegni and ManousakialJ find 
dW/t = 0.59 ±0.15 in the thermodynamic limit for d = 2 
and J/t = 5, while Eq. (§|) gives 0.41. 

Finally, we can extrapolate to two dimensions the char- 
acter of the ordered arrays of charged domain walls at 
low doping concentration and weak Coulomb interaction. 
Domain walls in two dimensions are one-dimensional 
(lines) and such ordered arrays are known as "stripe 
phases" . Directly extrapolating the optimal large dimen- 
sional domain- wall structures to d = 2, we would expect 



the stripes to be site-centered, vertical, antiphase domain 
walls in the antiferromagnetic order, and to be metallic 
(and possibly superconducting) for J/t < Y c and insu- 
lating for J/t > Y c . In particular, if we extrapolate the 
leading order expression for the electron density in the 
hole-rich phase, Eq. to d = 2, and then evaluate 

it for t ^> J, we find that such stripes should have ap- 
proximately 0.31 doped holes per site along the stripe, 
and are thus metallic. Transverse to the stripe direction, 
such a phase is a generalized charge and spin density wave 
state, in which the period of the charge density wave is 
half that of the spin density wave.Ej However, because 
of the electronic motion along the stripe, this phase is 
actually an electron smectic.0 Unfortunately, there are 
no detailed microscopic two dimensional calculations to 
compare with these results, so we compare them with 
experiments on doped antiferromagnetsEj 

C. Rigorous Results 

In addition to our perturbative results in powers of 
1/d, we have obtained rigorous upper and lower bounds 
on the ground-state energy of the undoped system. These 
bounds, which are also quoted in Table 1, are shown to 
converge in the limit d — ► oo. 

D. Relation to Experimental Results on Doped 
Antiferromagnets in Quasi-Two and Three 
Dimensions 

By now there are many examples of antiferromagnetic 
insulators that can be chemically doped. One promi- 
nent feature of these materials is the occurrence of high 
temperature superconductivity, a phenomenon ibr which 
the present results provide little direct insightO. How- 
ever various spin and charge ordered states, as well 
as "nearly ordered" fluctuating versions a£|,such struc- 
tures, have been observed in these systemsE3c3 by direct 
structural probes, especially neutron scattering. Two 
concrete, and well studied examples of this are the 
quasi two-dimensional perovskites La2- a; Sr : j : Ni04 and 
Lai.g-ajNdzSr^CuOj, in both of which the doped hole 
concentration is equal to the Sr concentration x. The 
undoped parent compounds (with x — 0) are antiferro- 
magnetic insulators with spin S = 1 for the nickelates 
and S = 1/2 for the cupraies. In both cases, upon dop- 
ing, the system formsESEj a "stripe" phase, in which 
the doped holes are concentrated in antiphase domain 
walls in the antiferromagnetic order. At present it is 
not known whether the domain walls are site or bond 
centered in general. (At higher doping concentration in 
the nickelates, there is strong evidence that both types 
of domain wall coexist due to interactions between the 
wallscj.) However, there is a crucial difference between 
the domain walls in the two materials: In the nickelates, 



4 



there is one doped hole per site along the domain walls, 
and the doped system is, correspondingly, insulating. In 
the cuprate, the hole concentration along the domain 
wall is roughly one doped hole per two sites along the 
domain wall, and the system is correspondingly metal- 
lic, and even superconducting, despite the presence of 
almost static charge and spin density wave order. (This 
latter behavioK-as very suggestive evidence of an electron- 
smectic phase!]) In addition, the domain walls are diag=. 
onal in the nickelatesEj and vertical in the cupratesB§EZl 
We feel that the occurrence of charged stripes in 
lightly-doped antiferromagnets, the fact that these 
stripes are antiphase domain walls in the antiferromag- 
netic order, and that they can be metallic or insulat- 
ing, depending on the ratio of J/t, are physically ro- 
bust features of the large d theory which we expect to 
apply mutatis mutandis in d = 2. However, the prefer- 
ence for vertical versus diagonal stripes, and site-centered 
versus bond-centered stripes is likely to depend on mi- 
croscopic details, even in large dimensions. Of more 
profound importance is the fact that, while in large di- 
mensions the charged domain walls always crystallize at 
low temperature into an ordered density wave, in low di- 
mensions, especially in two dimensions, there is the very 
real possibih.ty_±ha.t the domain walls will be quantum 
disordered .E3l3O'0 In such a melted state, which might 
be either fully disordered (isotropic) or still retain ori- 
entational order ("electron nematic"), the sort of charge 
and spin-ordered states that are characteristic of the large 
d theory occur as local correlations in the fluctuation 
spectrum; a microscopic electronic theory of such quan- 
tum disordered states is not available at present. 

II. THE MODEL 

The model we consider is the straightforward general- 
ization of the usual t — J model (or t — J — V modeltil) : 

<hj> 

<i,j>,a 

where Si = a , cj a o at(J iCi. a i is the spin of the electron 
on site i, rii = cj a c%,a is the number of an electron on 
site i, ct creates an electron with z-component of spin 
equal to a = ±1/2, a are the Pauli matrices, there is a 
constraint of no double-occupancy of any site, 

m = o, l, (2) 

and < i,j > signifies nearest-neighbor sites on the d- 
dimensional hypercubic lattice. In comparing results of 
different calculations, it is important to note that there 
is more than one definition of the t — J model. Most 



commonly,t2lli-3 "the t — J model" is defined as in Eq. 
(|]) with V = — J/4, but without the prefactor of 1/d. 
Where it can be done readily, we will quote results for 
arbitary V, but where this leads to complications, we 
will, for simplicity, analyze only the canonical case V = 
—J/4. The additional factor of 1/d is included so that the 
ground-state energy density remains finite in the d — > oo 
limit; thus, in making a comparison with previous results 
on the d — 2 t — J model, all energies computed here 
should be multiplied by d — 2. 

III. THE UNDOPED ANTIFERROMAGNET 

The undoped system has one electron per site so that 
the electron hopping term (t) has no effect, and the sys- 
tem is manifestly insulating; the only remaining degrees 
of freedom are described by a spin 1/2 Heisenberg anti- 
ferromagnet with exchange coupling J. 

A. Rigorous Bounds 

It is possible to obtain upper and lower bounds on 
the ground-state energy of the spin-1/2 Heisenberg model 
which approach each other in the large d limit. An upper 
bound is obtained by calculating the variational energy 
of the "Neel" state, which has alternating up and down 
spins on alternate sites, and gives a ground-state energy 
per site of En eel = —J/i + V. 

A lower, bound for the ground-state energy can be 
obtainedEII as follows: We express the full Hamiltonian 
as a sum of pieces, 

H= H 3 ( 3 ) 

j — black 

where the sum is over all sites on the black sublat- 
tice and Hj is the exchange interactions between site 
j and its nearest neighbors (which are necessarily on 
the red sublattice). The Hamiltonians Hj are read- 
ily diagonalizcd, but not simultaneously since they do 
not commute with each other. Nonetheless, the sum of 
the ground state energies of Hj gives the lower bound 
Eiower = — (1 + d^ 1 ) J/4 + V for the ground-state energy 
per site. 

These results, combined, prove that the ground-state 
energy per site, Eaf , of the Heisenberg model approaches 
that of the classical Neel state in the limit of infinite 
dimension, 

-(J/4)[l + l/d]<E AF -V <-(J/4). (4) 
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B. Perturbative expression for the ground-state 
energy and sublattice magnetization 

We now embark on the derivation of results in a sys- 
tematic expansion in powers of 1/d. For this purpose, 
we will consider the Heisenberg model as the isotropic 
limit of a Heisenberg-Ising model. To begin with, we use 
Rayleigh-Schrodinger perturbation theory to evaluate the 
properties of interest in powers of the XY coupling, and 
then reorganize this perturbation theory in powers of 1 jd. 
Thus, we take as our unperturbed Hamiltonian the Ising 
piece of the interaction, 



H = ^ ^2 [ J z s i s j + VriiUj] , 



<h3> 



and treat the XY piece 
Jl 



d 



E i s ? s l 



sis]] 



(5) 



(6) 



<i,3> 



as a perturbation. 

The ground state of H is the (two-fold degenerate) 
Neel state. Hi has the effect of flipping pairs of spins, 
which because of the large coordination in high dimen- 
sions means that the intermediate states have energies 
that are proportional to d. We have evaluated the per- 
turbative expression for the ground-state energy per site, 
Eaf, and the ground-state sublattice magnetization, m, 
to fourth order in y = Jj_/ J z , but it would be straightfor- 
ward (using modern methods of high temperature series 
expansion) to extend these results to higher order. The 
results are: 



m 



1 y y 

Ad U 2 
y 2 (128 + lV) 
1024d 4 



y 2 (^-y 2 ) 

256d 3 



+ 0(l/d 5 ) 



(10) 



The appropriate expressions for the Heisenberg model 
can now be obtained by taking the limit y — > 1. 



C. Goldstone Modes and the Long- Wavelength 
Physics 

Because the Neel state involves a broken continuous 
symmetry, we know that there must exist a gapless 
Goldstone mode, the magnon. In the presence of Ising 
anisotropy, the magnon is massive, and is perturbatively 
related to the single spin flip. Thus, we could imagine 
using the same decomposition of the Hamiltonian into 
an Ising and XY piece to compute the magnon spec- 
trum perturbatively, and then reanalyze the expression 
in terms of the 1/d expansion. This is impractical, but 
it is instructive to see why. 

In th order (i.e. in the Ising model), there is a set 
of N/2 degenerate excited states with excitation energy 
e = J z /2 and S z = — 1 obtained by flipping a spin on 
the black sublattice, and there is a complementary set 
of excited states with S z — +1 obtained by flipping a 
spin on the red sublattice. These states resolve them- 
selves into the two polarizations of the magnon band 
upon performing degenerate perturbation theory in pow- 
ers of y — J±/ J z . The results of degenerate perturbation 
theory can be summarized in terms of an effective Hamil- 
tonian, 



Eaf = V - 



1 



V 



y 4 (2rf-3) 



(2d 
+0(y 6 ) 



1) 4(2d- l) 3 (4d-3) 



and 



d[8d 3 + 2d - 3] 



^ r-i 2 ^ i 4 

171 = 2 [ ~ V (2d- l) 2 V 8(2d- l) 4 (4d-3) 2 (d- 1) 



+0(y 6 



(8) 



It is clear that each power of y brings with it an additional 
power of 1/d from the additional energy denominators, 
as promised, so that the 0(y 6 ) term is actually 0(y 6 /d 5 ). 
Reorganizing these expressions in powers of 1 /d yields 



Eaf 



1 y 2 y 2 y 2 (8-y 2 ) 

2d 4d 2 64d 3 



+ 



y 2 (16-3y 2 ) 
256d 4 



+ 0(l/d 5 ) 



(9) 



and 



H 



eff _ 



iSzj— black 



(ii) 



(7) 

where b\ creates a spin flip on site j and obeys boson 



ti 



commutation relations, [bi,b 3 



Sij . To be concrete, 



we have considered the magnon with S z = —1, so we 
take the Hamiltonian to operate in the 1 spin flip sec- 
tor, X^-^ocfc bjbj = 1. The effective Hamiltonian can 
be solved by Fourier transform to give a magnon energy 



mag(k) —J2j=black J ?-J ^ l"* ' 

interested in the case in which there was substantial Ising 
anisotropy, we could simply compute Ji.j to the desired 
order, since if i and j are n steps apart on the lattice, 
Ji,j ~ Jz[y/d] n , and hence for small y, H e ^ is short- 
ranged. It would also seem that the same logic would jus- 
tify the self-same expansion for large d, and indeed (as is 
implicit in the discussion of the ground-state energy) this 
is crudely true. However, even though J7ij falls rapidly 
with n, the number of n th "Manhattan" neighbors grows 
just as rapidly, i.e. as d n . For non-zero wave vector, this 
does not matter, as the far neighbors contribute to e mag 
with rapidly varying phases, and so the long-range tails 
of Jij arc unimportant. However, for k very near k = 



jToj exp [ik ■ Rj]. If we were actually 



G 



(or, equivalently, near k = tt =< tt, tt, tt, ... >), all terms 
in the Fourier transform add in phase, so Jij must be 
computed to infinite order. 

Of course the point is that the low-energy Goldstone 
modes have exceedingly small phase space in large di- 
mension, although they always dominate the temper- 
ature dependence of thermodynamic quantities at low 
enough temperature and the asymptotic decay of correla- 
tion functions at large enough distances. Thus the Gold- 
stone modes are entirely unimportant in high dimensions, 
except for physical quantities that strongly accentuate 
the lowest energy excitations. 

The way to study the Goldstone behavior is in terms 
of a spin-wave expansion, again suitably reinterpreted in 
terms of the 1/d expansion. We thus start by consider- 
ing the spin-S Heisenbecg. antiferromagnet in d dimen- 
sions using the standardo Holstein-Primakoff bosons to 
obtain the spin- wave spectrum in powers of 1/ S. We will 
confine ourselves, here, to the lowest order theory, as it 
adequately illustrates the point. The sublattice magne- 
tization is thus S in the classical Neel state, but receives 
a correction of order S° from spin-wave fluctuations as 

or, 1 f dk r 1 , 



+o{ Y2 )} 



(12) 



where the integral over k is over the first Brillouin zone, 

d 



7 



i%) = (dr^cos^j 



(13) 



a=l 



is the normalized structure factor, and the spin- wave en- 
ergy is 



e mag (k) = J5 V / [l- 7 2 (fc)][l + 0(1/ S)]. (14) 
Expanding the integrand in powers of 7 and employing 

dk 7 2 ™(£)/(27r) !i = (2n)!/[(4d)"ra!] (15) 



(16) 



we obtain 



Clearly, in the limit S = 1/2, the spin- wave expansion 
can be re-expressed as an expansion in powers of 1/d. 
Indeed this expression agrees with the earlier result of 
perturbation theory in y, in the limit y = 1, as, of course, 
it must. 

For k near (or equivalently, near if), the spin- wave 
spectrum (for S = 1/2) can be expanded in powers of |A;| 
to give the usual linear dispersion of the Goldstone mode 
with spin- wave velocity 



(J/d)y/d/2[l + 0(l/d)]. 



(17) 



We can also compute the transverse spin-spin correla- 
tion function to leading order in 1/S, and examine the 
resulting expression at large d. Using standard results, 
it is easy to seetS that 



< [sfs% + sfsy] >= se**-*** 

[l + 7(fc)]e"'^ + [I-7®] 



l-7 2 (fc) 



dk 

cos (k ■ Rij) 



[1 + 0%)]. 



(18) 



where Rij = R4 — Rj. In the large R limit, this integral 
can be evaluated by approximating the integrand by its 
small k expression, 1 — 7 2 (fc) ~ rf _1 /c 2 , which yields the 
Goldstone behavior, 



< [SfSf + S?S?] > 



5T(- 



2)V^f( 1 



25 



Ro 



R 



TT y/rTR 

(d-1) 



(19) 



where R = \Rij\, Ro — \J d/2~Ke, V is the gamma func- 
tion, and in the second line we have used Stirling's for- 
mula for large d. It is easy to see that the integral eval- 
uated above is dominated by values of k 2 <~ d 2 /R 2 , and 
since the small k approximation is valid only so long as 
k 2 <C d, the Goldstone behavior is only valid for R 3> Ro, 
as is suggested by the form of the result. 

The most efficient way to evaluate properties of the 
system at low but non-zero temperature is to use the 
1/d expansion to compute the fully-renormalized zero- 
temperature parameters that enter the 0(3) nonlinear 
sigma model which governs the Goldstone modes, namely 
the spin- wave stiffness, p s , and the transverse uniform 
susceptibility, \o- The susceptibility can readily be com- 
puted perturbatively in powers of 1/S, and the resulting 
expression reexpressed in powers of 1 / d as 



Xo 



1 

4X 



1 



8Sd) 



0( 



1 



d 2 S 2 



) 

(20) 

In terms of this, p s can be computed from the relationEl 
Ps = c 2 X o (21) 
where c is the spin- wave velocity given in Eq. ([Pfl). 



IV. ONE HOLE IN AN ANTIFERROMAGNET 

The one-hole problem has a structure that is nominally 
like that of the one magnon problem, but added factors 
of 1/d make the perturbative approach tractable for all 
values of k. We define as the unperturbed Hamiltonian 
Ho, the Ising limit of the t — J model, with J± = t = 0. 



7 



A. The minimal hole with S 2 = 1/2 

There are 2N degenerate one-hole ground states of Hq, 
where the factor of 2 is due to the global degeneracy of 
the Neel state, and the factor of N (which is the number 
of lattice sites) comes from the locations of the empty 
site (the hole). (Henceforth we focus only on the states 
in which the magnetization is up on the black sublattice.) 
These states are, in turn, separated into disjoint Hilbert 
spaces labeled by the conserved quantum number, the 
total z-component of spin, since a hole on a black sublat- 
tice site has S z = +1/2 and one on the red sublattice has 
S z = —1/2. For concreteness, we will focus on the N/2 
degenerate states corresponding to a hole on the black 
sublattice. 

We now use degenerate perturbation theory to con- 
struct the effective Hamiltonian of one hole, 



H 



i&cj — black i&zj—red 



where cj is the fermionic creation operator for a hole on 
site j, and for Hermiticity, — tji. Once the effective 
Hamiltonian is computed, its eigenstates and eigenvalues 
can be found by Fourier transform: 



ehoie(k) = 2J t j exp[ik ■ Rj]. 

j— black 



(23) 



To begin with, we study the perturbative expressions 
for the diagonal term, eo = tu ■ 



e = -2V- 



-2V 



Jz 



i + ir 



2(d-l) 



(2d-l)(4d-3) (2d-l) 



+0(y 4 ) + 0(z 4 ) + 0(y 2 z 2 ) 



(24) 



[y 2 - 16z 2 } [y 2 - 32z 2 ] ^ . - ' 
1 + + , „ — - + 0(l/d 3 ) 



16d 2 



where z = t/J z and y = Jj_/J z . The corrections to the 
hole self energy which are independent of y are the fa- 
mous string corrections to the hole self energy, of which 
the retraceable paths of Brinkman and RiceE3 are a sub- 
class. 

Next we study the coupling between second "Manhat- 
tan" neighbor sites (nearest neighbors on the black sub- 
lattice). There are 2d(d — 1) "true second neighbors", 
reached by taking a step to the nearest neighbor site in 
one direction, and then a second step in an orthogonal di- 
rection, and there are 2d "straight line" second neighbors 
reached by taking two steps in the same direction. For i 
and j true second neighbors, ty = 2r + rising, where the 
factor of 2 in the definition takes account of the fact that 
there are two minimal paths to the true second neighbor. 
Here r is given by 

ttyz 



d(2d- l)(4d-3) 



1-2/ 



(6d 3 - 6d 2 



,N 2 



D 2 
314d- 7) 



3(d- l)(2d-l)(4d-3) 
2 (528d 4 - 1368d 3 + 1160d 2 - 
~ Z 3(2d- l) 2 (d- l)(4d- 3)(6d- 5) 

+0(y 3 ) + 0(z 4 ) + 0(z 2 y) 



tyz 
2d 3 



(l 



4d' 



22z' 
^5" 



<D(l/d 2 



(25) 



with 



N 2 = 294912d 7 - 1834752d 6 + 5070688d 5 - 7924688d 4 

+7461296d 3 - 4195852d 2 + 1298748d - 170325 (26) 



and 



Do 



96(2d - l) 2 (d - l)(3d - 2)(4d - 3) 2 
x(6d-5)(8d-9)(8d-7); (27) 



r comes from processes in which, in lowest order, the hole 
hops twice, followed by a spin exchange which repairs the 
resulting damage to the spin order. In addition, there is 
a contribution 



32tz- 



Tl 



sing 



3d{2d- l) 2 (4d~3) 2 (d- 1) 



l + 0(z 2 ) 



tz° 
~6cF 



1 



(28) 



which comes from a process in which a hole circles a 
plaquette one and a half times, thus "eating its ow 
string" . This process, which was discovered by Trugma: 
for d = 2, survives even in the Ising limit y = 0. How- 
ever, Ti s i n g is higher order in powers of 1/d and so is 
negligible, even when t ^> J. For straight-line second 
neighbors, ty = r' + r' Ising where 



Atyz 



d(2d-l)(4d-3) 
(d-1) 



y 



(4d - 3) 



y 



2^2 

D> 



' D' 3 



0(y 3 

tyz 

2d 3 



(1 



0(z 4 ) + 0(z 2 y) 



4d 



(1 



22Z 2 
^5" 



0(l/d 2 ) 



(29) 



where 



N 2 = 73728d 6 



41088d 5 - 514544d 4 + 1145088d 3 
1016784d 2 + 421204d- 67611 (30) 



D' 



48(2d - l) 3 (3d - 2)(4d - 3) 3 (6d - 5) 
x(8d-9)(8d-7) (31) 



and 
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^3 



-176d 3 + 280d 2 -108d-2 



D' 3 = (2d- l) 3 (4d-3) 2 (6d-5) 



(32) 



(33) 



and Tj sing is the corresponding Trugman term which is 
0(tz 9 /d 9 ). Because they are higher order in 1/d, we will 
henceforth ignore Ti S i ng and T' Ising relative to r. How- 
ever, we shall see below that the difference, 



Atzy 2 (M- 2f 



3d(d- 
3tzy 2 
16d 4 



l)(2d- 
i , ( 13 



l) 2 (4d 



ad 



-3) 2 

0(dr 2 



l + 0(y)+0(z 2 ) 



(34) 



although smaller than r by a factor of 1 /d, plays a critical 
role in determining the band structure of the minimum 
energy hole in an antiferromagnet. 

The sign of these terms deserves some comment. Since 
the lattice structure defined by is not bipartite, the 
sign of the matrix elements is physically significant. In 
the present case, since the leading order contributions to 
r come from third order perturbation theory, the result- 
ing matrix element is positive. 

These expressions may be combined to obtain an ex- 
pansion for the bandwidth W of a single hole by adding 
the contributions from the different hopping processes. 
Each hop contributes a factor 2r or 2t' to W, so 



W = Ad(d - l)r + 4dr'. 



(35) 



Then, using Eqs. ( |25| ) and (|29j) and expanding in powers 
ofd-\ 

dW „ , y. 5yz ,„ y 22 z 2 . „ . , 



This result does not make sense unless W > or 



15 r Ad, 



<^r™(i_^) + i_^l 
22 L 5 V 4 7 2 J 



(36) 



(37) 



Since z = t/J z , this illustrates the fact that the motion 
of a single hole is exchange dominated in the large-d ex- 
pansion. 

It is important to note that the leading order expres- 
sion for r contains one more power of 1/d than the corre- 
sponding matrix element in the effective Hamiltonian for 
one magnon. Indeed, the leading-order behavior of the 
matrix elements connecting sites separated by 2n nearest- 
neighbor steps (2n th nearest Manhattan neighbors) is 



2)1 



t 2n ~ tz 2n - 1 y n /d 3n , where the factor of iz 2 ™" 1 oc t 
comes from the minimum number of hops for the hole to 
propagate this distance, the factor of y n oc J™ reflects 
the minimum number of spin flips needed to restore the 
spins to a ground-state configuration following the pas- 
sage of a hole, the factor of d~( 3n_1 ) comes from the 



accompanying energy denominators at this order of per- 
turbation theory, and one additional factor of 1 j d comes 
from the overall normalization of the Hamiltonian. Be- 
cause of these extra factors in the expression for , its 
contributions to the hole energy fall rapidly with distance 
in high dimensions, where the number of 2n-stcp paths 
is 2d(2d - l) 2n ~\ and so the number of 2n th Manhat- 
tan neighbors can grow no faster that (2d) 2n . Thus, the 
longer range pieces of ty can be neglected for any value 
of k in large enough dimension. 

The effective Hamiltonian obtained by retaining only 
terms out to second Manhattan neighbors is given by 

e h ole(k) = e - 2d(r' - r) + (2d) 2 T1 2 (k) 

d 

+4(r'-r)^cos 2 [fc a ]. (38) 

a=l 

If we ignore the small difference, (r' — t), then the 
minimum energy hole states are located on the d — 1 
dimensional hypersurface, j(k) = 0; since the band- 
structure for the noninteracting tight-binding model is 
e/ree = — ^tj(k), and that model is, in turn, particle- 
hole symmetric, this hypersurface is precisely the Fermi 
surface of the half-filled band in the absence of interac- 
tions. With higher order terms in powers of 1/d (which 
produce a non-zero value of (r' — r) > 0) the minimum 
energy of a single hole occurs at k = (1/2)7? and the 2 d — 1 
symmetry-related points. 

We emphasize that, although we have treated the ef- 
fects of t perturbatively, we have not made a small t 
approximation. The present results are valid for arbi- 
trary t/J, so long as it is not parametrically large (i.e. 
so long as t <gC dJ). Nonetheless, because each hop of the 
hole may break 0(d) bonds, the large-d limit is exchange 
dominated. 



B. Magnetic polarons with larger spin 

In low dimensions, and for t S> J, it is believed that 
a single hole in an antiferromagnet produces a ferromag- 
netic bubble in its vicinity, or, more precisely, that there 
is a series of level crossings as a function oft/J at which 
the total z component S z of the spin of the ground state 
for a single hole state steadily increases. However, in the 
large d limit, the antiferromagnetic energy always domi- 
nates unless t is parametrically larger than J, i.e. unless 
t ~ d x J, where x is a positive exponent (which we will 
estimate below). Such parametrically large values are 
beyond the scope of the present analysis. 

To estimate the magnitude of t/J at which ferro- 
magnetic bubbles first appear, we consider a straight- 
forward generalization (tO|-sf > 2) of the calculation of 
Emery, Kivelson, and Lir£3 for a hyperspherical ferro- 
magnetic polaron in the large size limit (where discrete 
lattice effects can be neglected). We balance the mag- 
netic energy lost in the volume of the polaron against 







the zero-point energy (computed in the effective mass 
approximation) to localize the hole in the interior of 
the polaron. This results in a polaron with a radius, 
L w [2tn/A d (V - E AF + J/4)] 1 /( rf + 2 ), where A d is the 
area of the unit d dimensional hypersphcrc, 



A d = 2(V^) d /r(d/2) « ^d/^(^2^/d) a 



(39) 



and in the final line, we have used Stirling's approxima- 
tion for large d; the polaron has spin S z w A d L d /d. By 
definition, the spin of the polaron must be substantially 
greater than 1, which means that t ^S> d 2 J. We conclude 
that in the large d limit, the low energy hole branch is 
always the naive, S z = 1/2 vacancy state, and that local 
ferromagnetism is never a relevant piece of the one hole 
physics. 



V. EFFECTIVE INTERACTIONS BETWEEN 
TWO HOLES 

Broadly speaking, the effective interactions between 
two holes are of two kinds - "potential" , which are in- 
duced by distortions in the antifcrromagnetic order, and 
"dynamic" , which minimize the zero-point kinetic energy 
of a hole. 

At long distances, the effective potential of interaction 
can be computed by considering change in the magnetic 
Hamiltonian induced by two static holes at lattice sites 
and i, 



Ho 



(i) 



Si 



(o) _ ^■ 

E So ■ Sk 



(40) 



where Y]j signifies the sum over the nearest neighbor 
sites of i. Then, on integrating out the magnetic degrees 
of freedom, we obtain 



W (o) 



V 



eff 



(Ri) = -J 2 /d [dtf^Y. 



3 k 



T[Si(t) ■ SjQb)- < Si ■ Sj > ] [So -S k -< S -S k > ] 
+ ... (41) 



where ... are higher order terms in powers of J, which are 
also of shorter range as a function of f is imaginary- 
time, and T is the-imaginary time ordering operator. It 
is straightforwardly to determine from linear spin-wave 
theory that 



(42) 



which is a short range potential in the sense that the inte- 
gral over space of |V e -^| is non- infinite in all dimensions 
greater than d = 1. (The integral over space of V e ^ 
itself is easily seen to be 0.) Moreover, the long distance 



tails of V e M decrease in importance as d increases. For 
this reason, we will ignore the power-law tails of V e ^ 
and simply consider its dominant, short-distance pieces. 

The nearest-neighbor interaction between two holes 
(which is certainly attractive in the canonical case, V — 
— J/4) can readily be computed from perturbation theory 
in powers of y = J±/J z : 



(43) 



There is a considerably weaker interaction (which, in fact, 
is repulsive) between second nearest-neighbor holes: 



V eff {e 1 +e 2 ) 



JzV 2 



Ad{2d 

J z y 2 



32d 3 



l)(4d-3) 
[l + 0(l/d)]- 



[l + 0(y 2 )] 



(44) 



Indeed, to this order, the effective interaction is the 
same for all second Manhattan neighbors, V e *' (2§i) = 
V efi {e x + e 2 )[l + 0(l/d 3 )}. Clearly, for further Man- 
hattan neighbors, the effective interactions are down by 
additional powers of 1/d. 

The kinetic terms, in general, generate fairly compli- 
cated interactions of the form 



T 



eff _ 



ijkl 



!; ic\c]c k ci 



(45) 



where, as before, c\ is a fermionic creation operator for 
a hole at site j. However, in large d, it is strongly domi- 
nated by its short-range components, of which the dom- 
inant terms are a potential interaction between nearest- 
neighbor holes (which renormalizes V e ^ (ei)) and a pair- 
hopping term. Indeed, combining the potential and ki- 
netic terms to leading order in 1/d, we find a two hole 
contribution to the effective Hamiltonian (i.e. the in- 
teraction part of the effective Hamiltonian, of which Eq. 
(P2|) is the noninteracting piece): 

H e 2 ff = u e " c hV* ( 46 ) 
<ij> 

~ T£ff E c]4c 3 c k + 0(l/d 3 ) 

<i,j,k> 

where in the pair-hopping term < i,j,k > signifies a set 
of sites such that i and k are both nearest-neighbors of j 
(which we define to include the case i = k), and 

A Art L 



d U [ 
+0(y i ) 



(d- l)(2d- 1) 
G(z 4 ) + 0(y 2 z 2 )} 



(47) 
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and 

= |[l + j(l + f) + 0(l/d 2 )]. (48) 

The pair- hopping term, T e ^ , has an interesting his- 
tory: In early work on high temperature superconductiv- 
ity, it was often claimed that, whereas the motion of a 
single hole is inhibited by antiferromagnetic order, pair 
motion appears to be entirely unfrustrated. It was sug- 
gested that this might indicate a novel (non-potential) 
source of an attraction between holes which could be the 
mechanism of high temperature superconductivity. At 
first sight, the fact that T e ^ ~ d~ 2 , while the single- 
particle hopping term is r ~ d~ 3 , appears to support 
the validity of this idea in large d. However, the fallacy 
of this argument was revealed in the work of TrugmanQ, 
who showed that this mode of propagation of the hole 
pair was frustrated by a quantum effect which originates 
in the fermionic character of the hole. In large d, this 
frustration effect is particularly graphic. Pair binding is 
enhanced if we ignore single-particle hopping, and diago- 
nalize H% . Of course, any state in which the two holes 
are farther than one lattice site apart are eigenstates of 
iff , so long as terms of this range are neglected (be- 
cause they are of higher order) as in Eq. ([47]) . For the 
states in which the two holes are nearest-neighbors, iff 
can be block diagonalized by Fourier transform, with the 
result that there are d bands of eigenstates labeled by a 
band index and a Bloch wave-vector, k. It is straight- 
forward to see that none of these bands disperses (their 
energies are independent of k) and that d — 1 of these 
bands have energy U eff , while the remaining band has 
energy U eff + 2T e H. This final band, which feels the 
effect of pair propagation, has the highest energy. On 
the other hand, if the holes were bosons, this latter band 
would have energy U e ^ — 2T e ^ , which is much closer 
to what one might have expected. 

It follows from this argument that coherent propa- 
gation of a pair is not an effective mechanism of pair 
binding and that the short-range attraction between two 
holes in an antiferromagnet arises from the fact that two, 
nearest-neighbor holes break one less antiferromagnetic 
bond than two far-separated holes. This interaction is 
sufficient to produce a two-hole bound state because the 
one-hole spectrum has an essentially degenerate band 
minimum along the d — 1 dimensional magnetic Brillouin 
zone. As in the Cooper problem, this gives a constant 
density of states at low energy and any attractive inter- 
action is sufficient to produce a bound state. 

VI. FINITE HOLE CONCENTRATION 

We have suggestedEBi that, in general, a doped anti- 
ferromagnet in d > 2 will phase separate into a hole-free 



antiferromagnetic region and a hole- rich regiop- .-There 
is now subs tant ial evidence, both numericalErE-j and 
analyticalEj'EZlH, that this is the case for the t — J model 
in d = 2. Phase separation is, of course, a first order 
transition, so it must be studied by comparing the total 
energy of various candidate homogeneous and inhomge- 
neous states to find the true ground state at fixed hole 
density. 

In large dimension, there are many metastable states 
which arc, in a sense, "local" ground states of given char- 
acter. While the large d limit allows us to compute the 
energy and character of a given candidate ground state 
exactly, it is almost never possible to prove that we have 
actually identified the global ground state. Specifically, 
we have computed the energy of various candidate states 
(as discussed below) and found that, of these, the lowest- 
energy state is phase separated into an undoped antifer- 
romagnet and a hole-rich phase with a very low electron 
density, i.e. with hole concentration equal to or nearly 
equal to 1. Moreover, given that the interaction between 
two holes is strongly attractive (of order the hole band- 
width) at short distances, and weakly attractive at long 
distances, we feel that it is extremely unlikely that any 
dilute hole-liquid or hole-crystal phase in the antiferro- 
magnet is stable in large d. Below, we show that the 
same instability shows up in a dilute domain-wall phase 
in which the holes are concentrated on an array of widely 
separated domain walls. 

What this means is that, under the assumption that 
we have not overlooked a lower-energy state (which we 
consider unlikely), we can obtain a complete and exact 
understanding of the zero-temperature phase diagram in 
the d —> oo limit by considering the phase coexistence 
between the undoped antiferromagnet (which we have 
already characterized) and a very hole rich or dilute elec- 
tron phase. 

We emphasise that the physics of phase separation and 
domain walls (Sec. VII) in large-d is not exchange domi- 
nated becuse it does not involve breaking a large number 
of bonds. 



A. Properties of dilute electrons in large d 

We now consider the ground-state properties of dilute 
electrons on a d-dimensional hypercubic lattice, for d 
large. Near the bottom of the band, the electron dis- 
persion relation is approximately quadratic, i.e. e(k) « 
— 2t + tk 2 /d is a good approximation so long as each 
component of k is small compared to 1, or typically, 
that k 2 <C d. Also in this limit, the interactions be- 
tween electrons are weak (since they rarely approach 
each other), and so can be ignored to first approxima- 
tion. Thus, we will begin by considering the properties 
of the non-interacting, quadratically dispersing electron 
gas in d dimensions. For this problem, the Fermi mo- 
mentum as a function of the chemical potential, [A, is 
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kp = */ (fi + 2t)d/t for [i > —2t, and the corresponding 
density is 



C. Effective Interactions in the Metallic State, and 
the Conditions for Superconductivity 



2A d fk F 
d V 2tt 



ird \ V 27rrf 



i-M + <>(?)] ( 49 ) 



where is given in Eq. ( p9| ) , a spin-degeneracy factor of 
2 has been included, and the final equality uses the large 
d expression for A^. Note that whenever the condition 
kp <C d is satisfied, the electron density is exponentially 
small for large d, so our approximations are exponen- 
tially accurate. The energy per site of this system can be 
computed readily: 



E, 



gas 



-2tn\\ 



\/(2d + A)]. 



(50) 



B. Conditions of thermodynamic equilibrium 

In general, for two phases to be in thermodynamic 
equilibrium, they must have equal chemical potentials, 
\x. However, here, the undoped antiferromagnetic phase 
is incompressible, so that the zero temperature chemical 
potential is undetermined. Then the condition for the 
electron gas to be in equilibrium with the antiferromag- 
net is 



fx = [E AF - E gas (ii)}/[1 - n(n)] 



(51) 



where E gas is the ground-state energy per site of the elec- 
tron gas and n is the electron density, both of which are 
funtions of /i. Thus, if Eaf < —2t, the only thermody- 
namically stable zero temperature phases of the t — J 
model are the undoped antiferromagnet and the vac- 
uum (no electrons). On the other hand, if Eaf > —2t, 
phase coexistence is possible between the antiferromag- 
net and a "metallic" phase with allowed electron densi- 
ties, n < n max , where n max is the electron density at the 
equilibrium value of /i. We shall see shortly that both 
E gas and n are exponentially small at large d, so that to 
exponential accuracy, 



fi = E AF . 



(52) 



If we use the large d expression for Eaf , we conclude that 
the metallic state is stable only if J < 4t[l + (D(l/d)], and 
that if this condition is satisfied, the maximum stable 
density of the metallic phase is 



e(4t - J) 



Airt 



+ (53) 



Notice that this quantity is small (and hence our approx- 
imations are justified), even in the limit t J, where 



(2/v^)[-] 



e i d/2 
7T J 



(54) 



Since the electron density in the metallic state is small, 
interaction effects are dominated by pairwise collisions 
between electrons. In the triplet channel, there is a 
nearest-neighbor electron-electron repulsion of strength 
(4V + J) J 'Ad, while in the singlet channel there is an infi- 
nite, on-site repulsion, and a nearest-neighbor attraction 
of strength {AV — iJ)/Ad. For the canonical choice of 
V = — J/4, which we adopt in most of this section for 
simplicity of notation, the attractive interaction in the 
singlet channel is simply —J/d. 

We can look for evidence of a simple s-wave instability 
of the metallic state in the low density limit by studying 
the conditions for the existence of a solution to the BCS 
gap equation. First, consider the unperturbed (noninter- 
acting) thermal Green function at low, but finite temper- 
ature 



G(k) = 



tanh[f(e(£)-/x)] 
2[e(£) - A 



(55) 



where e(k) — ~2t~/(k) and j(k) is defined in Eq. (|l3|). 
Now it is straightforward to show that, in the singlet 
channel, the BCS equation for the transition temperature 
T c may be written in terms of the corresponding real- 
space Green functions 



(56) 



a=l 



tN 



(57) 



and 



G2 = E g wE cos ( fc <*)] ; 



2t 2 N 



J2^ 2 (k)G(k) 



(58) 



as the U — > oo limit of 

= (1 + UG Q ){1 - JG 2 /d) + JUGf/2d 2 , (59) 
or, in other words, 

1 G 2 G\ 



J d 2d 2 G ' 



(60) 



Since e(fc) appears in the numerators of Gi and G2, as 
well as in the denominator of G{k), one may express G\ 
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and G2 in terms of Go and sums in which e(k) does not 
appear in the denominator. Now, to find the critical ratio 
of t/J at which T c — > 0, note that, for /x near the bottom 
of the band (/i rs —2t), the approximation tanh[/3(e(fc) — 
/i)/2] ss 1 can be used in the various sums, except in Go- 
Then Eq. ( J60| ) becomes 



J 2t 8i 2 G v ; 

Consequently the effective interaction is attractive and 
the BCS transition temperature for s-wave pairing is non- 
zero, so long as 

J/t > 2. (62) 

This result is valid for dilute electrons in any dimen- 
sion greater than 1, and is ia. agreement with earlier re- 
sults in two dimensions Notice that Eq. ([H]) with 
tanh[/3(e(fc) — fj)/2] set equal to 1 in Go is the condi- 
tion for a two-hole bound state, but in this case Go is 
divergent only in d = 1 and d = 2. Thus, in higher di- 
mensions, the critical value of J/t for a two-hole bound 
state depends on d, and tends to infinity as d — > 00. 

For J/t < 2, the direct pair interaction is repulsive in 
all channels, so there is no solution of the BCS equa- 
tions. However superconductivity could emerge in a 
higher angular-momentum channel by a variant of the 
Kohn-Luttinger effect. This sort of instability has been 
studied extensively in two dimensions ,Q and could prob- 
ably be analyzed in much the same way here. However, 
as the electron density is anyway exponentially small in 
large d, and these effects are of still higher order in the 
density, they occur at energy scales that are exponen- 
tially smaller than the exponentially small Fermi energy. 

One can, of course carry through the same analysis for 
arbitrary values of V , in which case, Eq. ( |si|) is replaced 
by 

3 J 

— -V>2t (63) 

and, in general, there is no superconducting transition 
if V > 3 J/4 - 2t. Recently, Riera and DagottcO have 
shown that a sufficiently large value of V prevents bound 
states of pairs of holes in the two-dimensional t — J model. 

Because there are no particularly good nesting vectors, 
we consider it unlikely that the hole-rich phase is subject 
to any sort of charge or spin-density wave instability of 
the hyperspherical Fermi surface in large dimension. 

VII. ONE DOMAIN WALL 

As mentioned previously, one feature of this large-d 
perturbation theory is that there are many interesting 
states which, if not the global ground state of the system, 
are at least the lowest-energy eigenstates in a restricted 
sector of the Hilbert space. In particlar, domain-wall 



states are interesting in their own right, although they 
prove to be of special importance in the presence of long- 
range Coulomb interactions. 

First of all, we define a domain wall to be a d — 1 di- 
mensional hypersurface which cuts the full lattice in two, 
such that far from the domain wall the system is un- 
doped, and hence antiferromagnetically ordered. Since, 
as we shall see, in large d such domain walls are always 
smooth (transverse quantum fluctuations of the position 
of the domain wall are suppressed by powers of 1/d) , we 
can characterize such a state by the mean position of the 
domain wall, its net charge (i.e. the concentration of 
holes per hyperarea) and the phase shift (if any) suffered 
by the antifcrromagnetic order in crossing the domain 
wall. In addition, since domain walls lie along lattice 
symmetry directions in all cases of interest, we can clas- 
sify them by the broken crystal point group symmetry. 
Here we analyze what we believe to be the physically 
most important domain walls, but we have not yet car- 
ried out an exhaustive study of other possible domain 
wall structures. 



A. "Vertical", site-centered domain wall 

A vertical domain wall lies perpendicular to a princi- 
pal axis of the hypercube, which we call the "x-axis" , and 
its location is thus specified by a single x coordinate. If 
the domain wall preserves reflection symmetry relative 
to a hyperplane that is perpendicular to the x axis and 
crosses the x-axis at a lattice site (which we will take to 
be x = 0), the domain wall is said to be site-centered; 
if this reflection plane passes half-way between two ad- 
jacent lattice sites (which we will take to be x = and 
x = 1), the domain wall is said to be bond-centered. We 
begin with the case of a vertical, site-centered domain 
wall with charge density of one hole per site, i.e. one 
that corresponds to a hypersurface of empty sites. 

We thus consider as the unperturbed part of the Hamil- 
tonian, Hq the Ising part of the t — J model with t = 0, 
but possibly with different choices of the Ising z-axis to 
the right (x > 0) and to the left (x < 0) of the domain 
wall. Then H± is all the remaining interactions, including 
all terms involving t. 

The ground state of Hq in the appropriate charge sec- 
tor is four-fold degenerate: it has no electrons in the 
domain wall (one hole per site) and one of the two possi- 
ble Neel states in each of the two "halves" of the system. 
The energy per domain-wall site (relative to the energy of 
the uniform, undoped Neel state) can now be calculated 
perturbatively, and it can be seen straightforwardly that 
the perturbation theory is again controlled by large d; for 
instance, from second order perturbation theory we find 

_ (d+l)J z f y 2 (4d 2 -3d-2) \ 

Cvert ~ 2d 1 + (d+l)(2d-l)(4d-3)J 
8J z z 2 

~ [2d - 1 + cos(S)] + "' 
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where S is the phase shift of the Neel order across the 
domain wall. 

It is clear that to this order in l/d, the energy is inde- 
pendent of 5. However, from the perturbative expression, 
we can extract the leading order ^-dependence: 

e ver t(5) - e vert {-K) = |^ [C0S( ^ + 1] |l + 0(l/d)\. 

(64) 

Clearly, then, the energy is minimized by an "antiphase" 
domain wall, with S = tt. The preference for an an- 
tiphase domain wall is produced by the transverse fluc- 
tuations of the domain wall, and it is likely to survive in 
all lower dimensions. The only cost of large transverse 
amplitude fluctuations of an antiphase domain wall is an 
increase in the surface energy, whereas, for any other 5, 
there are regions of impaired spin correlations. Indeed, 
we know from the solution of the one-dimensional elec- 
tron gas with repulsive interactions that, near half filling, 
the holes are solitons that act locally like antiphase do- 
main walls in that the magnetic correlations are shifted 
by 7r as one passes a hole, although of course there is no 
long range magnetic order in this case. 

The empty domain wall we have been considering is lo- 
cally stable so long as J/t is sufficiently large (J/t > Y c = 
4[1 + 0(1/ d)] ), but is unstable for smaller values of J/t. 
To see this, consider the state in which one electron is 
transferred from the surrounding antiferromagnet to the 
bottom of the domain-wall band. In first order degener- 
ate perturbation theory in t, the bottom of the domain 
wall band is easily seen to be e = — 2t(d— l)/d+ i.e. 
equal to the band bottom of the dilute electron phase to 
leading order in l/d. This means that, for the domain 
wall to be in local equilibrium with the surrounding an- 
tiferromagnet, it must have approximately (up to higher 
order corrections in l/d) the same electron concentra- 
tion, n max in Eq. ( |53| ) , as the hole-rich phase which can 
exist in equilibrium with the antiferromagnet. It is in 
this sense that a domain wall can be viewed as a form of 
local phase separation. Since even when J/t << 1, n max 
is exponentially small, the existence of dilute electrons 
within the domain wall when J/t < Y c does not sig- 
nificantly affect any of the other calculations described 
above. However it does imply that these electrons render 
the domain wall metallic and, under appropriate circum- 
stances, superconducting. This, of course, implies a qual- 
itative difference in the electronic properties of domain 
walls for small and large J/t. 

B. Bond-centered vertical domain wall 

If we fix the hole density at one (or approximately one) 
hole per hypersurface unit cell, as above, it is immedi- 



ately clear to th order in l/d that a bond-centered do- 
main wall will have considerably higher surface tension 
than a site-centered wall: On the one hand, more antifer- 
romagnetic bonds are broken by the bond-centered wall. 
On the other hand, only a very small fraction of the states 
in the electronic band have energy of order —td°, so even 
if we ignored the constraint of no double occupancy, and 
allowed the remaining electrons in the domain wall to 
minimize the kinetic energy part of the Hamiltonian, the 
gain in kinetic energy will be higher order in l/d than 
the loss of exchange energy. 

Of course, if we were to roughly double the concen- 
tration of holes, then a bond-centered domain wall can 
be viewed as simply two nearest-neighbor site-centered 
domain walls. This situation will be considered, below, 
when we consider interactions between domain walls. 



C. Other Domain Walls, Domain-wall Kinks, etc. 

There are many other kinds of domain walls. For in- 
stance, one could consider a "diagonal" domain wall, ei- 
ther bond- or site-centered, which is infinitely extended 
in d — 2 directions, and lies along a 45° angle relative 
to the two remaining lattice directions, which we will 
call "x" and "y" . As an example, one could construct 
a bond-centered diagonal domain wall by placing holes 
in the x-y plane along a "staircase" of nearest-neighbor 
sites obtained by first taking a step in the x direction, 
then a step in the y direction, and so on. To zeroth or- 
der, this diagonal stripe has the same energy per hole as 
the vertical, site-centered stripe. 

To compare the energies of different sorts of domain 
walls, we imagine finding the state of minimum energy of 
a system of size 2N x 2N in the x — y plane, and of infinite 
extent in the remaining d — 2 directions. Considering 
the projection of the problem on the x — y plane, we 
study the possible ground-states of the system with 2N 
holes per plane, in the presence of a strong staggered 
field on the boundary that favors an up spin on the red 
sublattice on the 2N sites nearest the lower corner (i.e. 
for points x = 0,0 < y < N and < x < N, y = 0) 
and the opposite field, which favors up-spins on the black 
sublattice, on the AN— 2 sites which form an upper "cap" 
i.e. the points x = 0, N < y < 2N, < x < 2N,y = 
2N - 1, and x = 2N-l,N <y < 2N. These boundary 
conditions force the system to have an antiphase domain 
wall of length 2N sites (or greater), but permits it to 
choose whether to have a vertical, diagonal, or a piecewise 
vertical domain wall with some concentration of right- 
angle kinks. 

The kink energy is readily computed perturbatively in 
powers of y and z, and it also can be re-expressed in 
powers of d : 

_ J z (y 2 (4d 3 -13d 2 + 12d-2) 
E H nk - -J\^ d _ 1 ^2d-3)(4d~5) 
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Az 2 (d-'S) 



d(2d-3)(d- 1) 



2d 



(y 2 + 6Az 2 ) 
Id 2 



+ 0(^)\, (65) 



which is manifestly positive in large d. (Here E^^ is the 
energy per site of the (d — 2) dimensional hyperline at 
which a site-centered vertical domain wall makes a right- 
angle bend.) Similarly, the energy of a zig-zag diagonal 
domain wall can be compared to the energy of the vertical 
wall computed above: 



^diag ^vert 



Md+1) 

Ad 
16z 2 



dy 2 



(d 2 - l)(2d- 1) 



{Ad - 3)(2d - 3) (2d - l)(d 2 - 1) 
+ 0(z 2 y ) + 0(y 4 ) + 0(z 4 ) 



Jz_ 

Ad 3 



(66) 



Thus it seems that the vertical domain wall has the 
lowest energy in large dimension. However, it is easy to 
see that this is a model-dependent result. For example, 
one can readily construct models on a "Cu-O" latticeo 
rather than a hypercubic lattice in which a diagonal, or 
bond-centered vertical domain wall has the lower energy. 



D. Interactions between two domain walls 

Two domain walls attract each other at long distances 
through the exchange of spin waves, in much the same 
way as two static holes. (See above.) Again, in high di- 
mensions, we expect this effect to be less important than 
the short-range attraction between domain walls. For in- 
stance, to leading order, there is an attractive interaction 
energy per unit hyperarea between nearest-neighbor ver- 
tical site centered domain walls which to leading and next 
to leading order is simply equal to the nearest-neighbor 
attraction between two holes, Eq. (E3|) above. 



VIII. 



BEHAVIOR IN LARGE BUT FINITE 
DIMENSION 



We have found that, in some instances, the electron 
kinetic energy t may play a relatively small role in the 
physics in large d, because only states exponentially near 
the band minimum have energies of order —tdr, while 
the bulk of the states have energies of order t/\[d. Thus, 
these states will only come into play when t/J gets to 
be parametrically large: t/J ~ Vd, where the large d 
theory is more difficult to control. In such a regime our 
results are less complete, and more subject to worries 
that there could be states we have missed. For instance, 
the perturbative treatment of the one-hole problem is no 
longer well controlled by large d: As discussed above, the 



effective hopping matrix element to the 2n th Manhattan 
neighbor is of order t(t / J) 2 ™ -1 / d 3n , while the number of 
such neighbors increases as d 2n . Thus, for t/J ~ yd, the 
contribution of far neighbor hops to the hole energy for 
k near or 7? does not decrease with n. Since this only 
matters for a very small fraction of one hole-states, while 
for generic values of fc, only the small n terms are impor- 
tant, it is unlikely that this problem leads to any signifi- 
cant changes in the qualitative physics of the dilute-hole 
problem. It does, however, mean that we cannot be quite 
as confident of the completeness of our understanding of 
the problem in this limit, as when t/J is not parametri- 
cally large. 

Nonetheless, with certain plausible assumptions and 
some guidance from the results of various studies in d = 
2, we can elucidate much of the behavior of the t — J 
model in this region of parameters, as well. (Note, as 
mentioned previously, our results here are consistent with 
those obtained using a somewhat different approach, for 
the large d Hubbard model. □) 



A. The Phase Diagram for Jd 1/2 /t > 1 

To begin with, we consider the behavior of the system 
when J/t ~ 1/Vd is parametrically small, but J\fd/t >> 
1 is still large. 

For energies away from the band edge, i.e. for e ~ 
t/Vd 7 the density of states per spin polarization of the 
tight-binding model on the hypercubic lattice in large 
dimension is readily computed by the method of steepest 
descents: 



p(e) = (27T)- 1 J dxe lxe [j {2tx/d)] 

= ^/fexp[- £ 2 d/4t 2 ] \l + 0(e 2 ) 



(67) 



From this it is clear that, so long as the electron den- 
sity is low enough, we can approximately ignore interac- 
tions between electrons, (i.e. the infinite, on-site U and 
the nearest-neighbor antiferromagnetic interactions), the 
density as a function of chemical potential (for fi < 0) is 
readily seen to be 



n(/x) = erfc(Vrf/i 2 /4t 2 ) + <D(n(fi) 2 ), 



(68) 



where erfc is the complementary error function, and the 
energy density is 



2t 



cxp[-dn 2 /At 2 }+0(n(n) 2 ). (69) 



Thus, there is a regime of parameters, 1 ^> J/t ^S> 
2/\fd in which the density of electrons is small (but not 
exponentially small) and interactions between electrons 
in the "gas" phase can still be neglected in any total en- 
ergy calculation. Under these conditions, \E gas \ <C |£^4f| 
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and 1 — n w 1, so from Eq. (pl|), it is easy to see that 
the density of electrons in a hole-rich phase in thermo- 
dynamic equilibrium with the undoped antiferromagnet 
is 

n max = erfc(jVd/4t) + 0(n 2 max ). (70) 



B. The Phase Diagram for 1 > Jd 1/2 /t 

When J/t is reduced still further, so that Jyd/t gets 
small, n max approaches 1, and it is no longer possible to 
ignore the effects of interactions in the hole-rich phase. 
In this limit, we lose the possiblility of quantitatively re- 
liable results based on our large-c? approach. However, 
for very small J/t and densities near 1, it is reasonable 
to expect the "electron gas" state to be ferromagnetic, 
at least locally. The ferromagnetic phase is noninteract- 
ing in any dimension when V = —J/4, corresponding 
to the canonical definition of the t — J model. In that 
case, the equilibrium between the ferromagnetic hole-rich 
phase and the undoped antiferromagnet can be computed 
exactly. Using the large d expression for the density of 
states in Eq. (|67]), we obtain the implicit expression for 

Mi 

fi=[E A F-Ef erro ]/[l-n(iJ,)], (71) 

where the ground-state energy of the ferromagnetic 
Fermi-gas is 

Eferro = 1= CX P [- d^/At 2 ] , (72) 

V 7T« 

the electron density is 

= (1/2)[2G( M ) - sign(/i)erfc(V^74t 2 )], (73) 

and 9 is the Heaviside function. We can evaluate this 
expression in the limit of small J\fd/t: 

2t I 

p = — J ln[2t/VndJ[l + O(jVdft)} (74) 
Vd 

and 

1 - n max = ^~[l + O(jVdft)} . (75) 

This, finally, corrects the unphysical aspect of the 
phase diagram for d — > oo, in that it implies that the 
boundary of the two-phase region approaches the zero 
doping axis (as {Jy/d/t)(ln[2t/V^dJ]y 1/2 ) as J/t -> 0. 
This is shown in the phase diagram in Fig. 2. (Note that, 
because all interaction effects vanish in the ferromagnetic 
state, the location of the coexistence line between the fer- 
romagnetic Fermi liquid and the undoped antiferromag- 
net can be computed accurately in any finite dimension 
d, for which E af i s known.Ej For d = 2, for instance, 
1 — n max — B\J J/t where B ss 0.61.) 

For parametrically small J/t and larger electron con- 
centration, the nature of the phase diagram in large but 
finite dimension is currently unexplored. 



IX. THE EFFECT OF COULOMB 
INTERACTIONS 

The t — J model has been widely studied because 
it is supposed to represent the most important low- 
energy physics of a system of strongly-interacting charged 
particles. It is assumed that the long-range part of 
the Coulomb interaction can be ignored provided it is 
fairly heavily screened by a surrounding dielectric back- 
ground. But this assumption is not valid for a state which 
is macroscopically inhomogeneous. In the presence of 
Coulomb interactions, we need to do thermodynamics at 
fixed mean particle density, and the system must be neu- 
tral at long length scales, i.e. phase separation is forbid- 
den. In a system with an average electron concentration 
1 < n < n max , and a long-range but "weak" Coulomb 
interaction in addition to the strong short-range interac- 
tions of the t — J model, we. encounter a class of phe- 
nomena that we have namedtl "frustrated phase separa- 
tion" . Here, the system is homogeneous (neutral) on long 
length scales, but inhomogeneous on short length scales, 
with interleaving regions that look locally like the two 
phases that would coexist in the absence of the Coulomb 
interactions. It is the purpose of the present section to 
explore the consequences of frustrated phase separation 
in the t — J model plus "weak" long-ranged Coulomb 
interactions in the limit of large d. 

We define the Coulomb interaction in d dimensions by 
a generalized Poisson equation 

-V 2 (b(r) = V-E=(Q/d)p(r), (76) 

where <fi(f) is the scalar potential, E(f) is the electric 
field, p{r) is the particle density, 

U(r) = -{Q/2d)p(r)^{r) = dE ■ E/2Q (77) 

is the energy density, and Q is the effective charge (or 
background dielectric constant) which determines the 
strength of the Coulomb interaction. One could, of 
course, imagine different ways of scaling Q in the large d 
limit. Given that we are after the physics of frustrated 
phase separation, we wish to take the limit in such a 
way that i) macroscopic phase separation is forbidden 
but ii) for a homogeneous state, the long-range part of 
the Coulomb interaction is unimportant compared to t 
and J. This is accomplished by taking the limit in such 
a way that Q does not depend on d. 

A. An effective model for "frustrated phase 
separation" 

In a previous publication, we considered a two dimen- 
sional model which we argued represented the physics 
of frustrated phase separation at a coarse-grained level. 
This model is easily generalized to arbitrary dimension: 
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B = ~ X>i-l)fo-l) 

<hj> 

0-/2)^Vc(i,j)[ai-dr][jj-dr] 



(78) 



where Oi = 1 if "site" i is a hole-free region and <Tj = — 1 
if "site" j is a hole-rich region, J' is a short-ranged "fer- 
romagnetic" interaction, which promotes macroscopic 
phase separation of the two coexisting phases, Vc is the 
Coulomb interaction, suitably defined on the lattice, and 



TV' 



(79) 



is the mean charge per "site." In this model, we imagined 
that sites represented small regions which were nonethe- 
less large enough that the local state could be described 
as being one of the two phases that would be in equi- 
librium with each other in the absence of the Coulomb 
interaction. In the present large dimensional context, it is 
possible to derive this effective Hamiltonian microscop- 
ically, identifying the sites in the effective model with 
the original sites in the t — J model, the Ci = 1 state 
with a site occupied by an electron, the cti = — 1 state 
with an unoccupied site, and J equal to the nearest- 
neighbor attraction between two holes, derived in Eq. 
( [43| ) above. This model is insensitive to the fact that 
the hole-rich phase has a non-zero electron concentration 
for J/t < Y c , but since the electron concentration is al- 
ways exponentially small, this error makes no difference 
in the energetics and structure of the various phases of 
frustrated-separation. Similarly, it ignores the fact that 
in each disconnected region of (holc-frcc) antiferromag- 
net, there is a potential ground-state degeneracy asso- 
ciated with spin-rotational symmetry; this might affect 
the finite temperature behavior of the system, but has 
no effect on the ground-state phase diagram. 

We have studiedcll the ground-state phase diagram of 
this model for d = 2 (and for V c (r) = Q/r, rather 
than the true two-dimensional "Coulomb" interaction, 
Vc{r) — Q l°gM)j-t|he results there, as well as the results 
of exact solutionis of a large N, "spherical" , version of 
the same model, lead to the conclusion that in all di- 
mensions, there are a number of ubiquitous characteris- 
tics of the the ground-state phase diagram of Eq. (78). 
Specifically, for very large Q/J, the ground state is the 
"Wigner crystal" phase, or in other words, the ground- 
state of the Coulomb term, itself, which, for a < is a 
fully (i-dimensional crystal of isolated sites with a = — 1, 
in a sea of sites with a = +1; for example, for a = 0, 
the "Wigner crystal" is a state in which one sublattice 
is occupied by a = +1 and the other by a = — 1. Con- 
versely, for Q/J very small, the ground states consist 
of a sequence of "stripe" phases of varying period (as 
a function of Q/J'), or in other words phases in which 
the density of a = —1 is a function of one coordinate, 
and independent of the remaining (d — 1) coordinates. In 
this regime of the phase diagram, for fixed a, the smaller 



Q/J, the longer the period, as discussed below. Between 
the striped phases at small Q/J and the Wigner crys- 
tal phase at large Q/J, there typically occur a sequence 
of more complicated phases that interpolate between the 
two extremes. In two dimensions, we found that these 
phases occupy an exceedingly narrow sliver of the phase 
diagram, but we do not know how generic this behavior 
is. 



B. The properties of the stripe phases 

If we confine ourselves to considerations of stripe 
phases, then the Hamiltonian in Eq. ([78]) can be reduced 
to an effective one dimensional model by summing over 
the values of the Ising spins in the (d — 1) dimensions 
perpendicular to the modulation direction: 

H stnpe = -{J/±d)Y^{Pj ~ - 1) (80) 

3 

+(Q/2d)J2\i-j\[a i -a][a j -a}. 

i,3 

The Madelung energies in this equation involve infi- 
nite sums, which are readily carried out numericallyuil, 
but cannot be done analytically. However, we can make 
rather good estimates by replacing the lattice sums by 
integrals. Specifically, for fixed hole concentration, x — 
(1/2)(1 — a), an array with period L, which consists of 
alternating stripes of <r = — 1 of width W = (L/2)(l — a) 
and intervening regions of width [L — W) of a = +1, is 
seen in this way to have energy per unit volume (defined 
to be volume associated with a single lattice site) 



E 



stripe 



-JWx 



QL 2 x 2 {l-xf 
24 



(81) 



This expression is readily minimized with respect to L, 
or equivalently the stripe-width W , with the result 



W = 



12J 



Q {l-xf 



(82) 



Finally, recalling that because the lattice is actually dis- 
crete, the allowed values of W are actually nearest integer 
approximants to this expression, we obtain the approxi- 
mate condition that width n > 1 stripes are stable for 



(2n + l) > 



2AJ 



Q (i 



> (2n-l), 



and that width 1 stripes are stable for 
8J x 



1 > 



q (i-xy 



(83) 



(84) 



Of course, for large enough Q/J, the width 1 stripe phase 
should give way to the other phases, mentioned above, 
and eventually at very large values to the Wigner Crystal. 
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An interesting aspect of this model is that, although 
it only explicitly involves the enumeration of the charge 
ordering (i.e. which sites are occupied and which are 
unoccupied), given the fact that there is an energetic 
preference for antiphase ordering of the spins across a 
charged stripe, we can reconstruct the ground-state spin 
order (up to a global rotation) by imposing the constraint 
that within regions of occupied sites (a — 1), the spins 
are antiferromagnetically ordered and, across unoccupied 
sites (a = —1), the antiferromagnetic order suffers a tt 
phase shift. (At very small Qj J , where the width of the 
charged stripe becomes greater than 1, a simple general- 
ization of the above microscopic calculations shows that, 
in large d, the correct ground-state energy and spin order 
may be obtained by viewing a thicker stripe as a collec- 
tion of nearest-neighbor fundamental (width 1) stripes 
and assigning a it phase slip per stripe.) 



X. GENERALIZED SPIN LADDERS 

There has recently been interest in the properties of 
spin- 1/2 Heisenberg "ladders" in d = 2, where a ladder 
is effectively a one-dimensional system which has finite 
width in all directions, save one, in which it is infinite. 
We can readily apply our analysis to the large d gener- 
alization of these ladders. For instance, we consider a 
generalized "two-leg" ladder, in which the lattice has a 
width 2 in d — 1 directions, and is infinite in 1 direction. 

Proceeding as above, we first obtain a rigorous upper 
bound, Eisi ee i = —(1 + l/d)J/8, and a rigorous lower 
bound, Ei ower = —(1 + 3/d)J/8 on the ground-state en- 
ergy per site. It is interesting, in this context, that the 
ground-state energy in the large d limit approaches that 
of the classical Neel state, even though this is a one di- 
mensional system, so we know rigorously that there is 
no true long-range magnetic order. Indeed, following the 
Haldane conjecture,!^ it is clear that for any finite d, this 
system will have exponentially falling magnetic correla- 
tions and a spin gap. However, this physics will only be 
manifest at very long distances (probably exponentially 
long in the large d limit), and at shorter distances the 
system will appear ordered. 

Again, without repeating the earlier analysis, we can 
compute the ground-state energy in perturbation theory 
in powers of y, which, to second order, gives 



Eaf = 



(d+l)J z 
8d 



1 + ^t 



y i (d 2 + 3d - 2) 
4d 3 (d + l)(2d- 1) 



0(y 6 ) 
(85) 



from which we can deduce that 



Eaf = ~ 



1 



J_ 

8 

16d 4 



(i + z/ 2 ) 

d 

0(l/d 5 ) 



d 2 



y 

8d 3 



XI. EXTRAPOLATION TO D = 2 AND D = 3: 

While the 1/d expansion gives us a small parameter 
with which to analyze the problem, it is legitimate to 
question the relevance of the 1 jd results for the physically 
interesting dimensions d = 2 and d — 3. (As there exist 
good methods for solving the present class of problems 
in d = 1, we are not concerned with pushing our results 
all the way down to d = 1.) 

For most of the types of order that we have considered, 
d = 1 is the lower critical dimension (for quantum disor- 
dering), and it is reasonable to expect large d results to 
be qualitatively reliable for zero temperature properties 
of the system in d = 2 and 3, but not for finite temper- 
ature properties in d — 2. As mentioned previously, this 
expectation is borne out to a large extent by comparison 
of the large d phase diagram of the t — J model shown 
in Figs. 1 and 2, with the best available data based on 
analytic and numerical studies of the model in d — 2. 

We can get more ambitious, and consider how good 
the quantitative agreement is between available numer- 
ical and series results for the model in physical dimen- 
sions and the results of the large dimension expansion. 
In Tables 1 and 2 we compare known results for the un- 
doped, spin 1/2 Heisenberg model in physical dimensions 
with the results of straightforward perturbation theory 
(in powers of y = Jj_/J z ) and of the 1/d expansion. 
Clearly, both give quantitatively excellent results. How- 
ever, whereas perturbation theory gives its best results if 
terms only to second order are retained, the 1/d expan- 
sion appears to approach closer to the correct value with 
each successive order, at least to fourth order (which is 
the highest order we have computed). In a sense, for an 
asymptotic expansion, it is the question of to how high an 
order do the results improve, even more than the overall 
accuracy of the result, which addresses the issue of how 
small is the expansion parameter. 





Eaf (2) 


m(2) 


Eaf (3) 


m(3) 


E2-leg{2) 




-0.5 


0.5 


-0.75 


0.5 


-0.375 


y 2 


-0.6667 


0.3889 


-0.9 


0.44 


-0.5625 




-0.6657 


0.3929 


-0.8995 


0.4404 


-0.5729 


Exact 


-0.6690 


0.3070 






-0.57800 



Table 2 : Comparison of the results of exact numerical 
studies (the row labeled "exact" ) on the two-dimensional 
spin-1/2 Heisenberg antiferromagnet, with the perturba- 
tive results in powers of y = J±/J z derived in the present 
paper. The dimension is indicated by the arguments of 
the computed quantities. The approximate results are 
obtained by setting y = 1, V = 0, and d = 2 or 3 in the 
series expansion, evaluated to the stated order. All ener- 
gies are measured in units of J/d, and the magnetization 
m is quoted in units in which g/is = 1, where /i_b is the 



(86) 
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Bohr magneton. 
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